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S '■ Abstract 

J>-»| Using experimental longitudinal and transverse velocities data for very high Reynolds number 

turbulence, we study both anisotropy and asymmetry of turbulence. These both seem to be 
. related to small scale turbulent structures, and to intermittency. We may assume that the large 

scale velocity shear gives an impact into the small scale turbulence, resulting in non-locality, and 
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INTRODUCTION 



Recently considerable attention has been devoted to the study of local isotropy of the 
high Reynolds number turbulence suggested by Ij. The turbulence is stirred at the large 
scales, and this energy is transported into the small scales, where, after intense nonlinear 
interaction, it becomes isotropic. Deviation from isotropy would mean that there is direct 
interaction between large scales containing non-universal anisotropy and small scales, leading 
to non-universal behavior of small scale spectral properties. 

It was indeed shown experimentally that in a sheared turbulence the isotropy is not 
.fficently .«.o.d fo. both .ca,a. and veloCy Eeld. fl. g g g g g J 

High Reynolds number (atmospheric) turbulence also shows deviations from isotropy. 
It was shown that the large scale shear does contribute into the scaling of the structure 
functions, see Anisotropic scaling of high-order structure functions was studied by 10 1. 



It was shown that the anisotropy in small scales remains stronger than ex pec ted before. 



The SO (3) decomposition was used to describe the anisotropy, see [lO|], [HI], [121], [iSj. The 
persistent anisotropy in small scale turbulence was found to be related to the intermittency 
corrections, 

In numerical simulations, the failure to return to isotropy was linked to both asymmetry 



of the probability distribution function (PDF) and to the vortex sheets, 15j. It became 
clear that the shear in the integral scale induces asymmetry down to the small scales, where 
it is manifested by intermittent structures like cliffs, etc., |l6| . 

On the other hand, the asymmetry PDF was found to be related to the intermittency. 
Thus, these three items, i.e., anisotropy, asymmetry and intermittency seem to be 
related. Note, however, that, in principle, these items are independent of each other. For 
example, the asymmetry of the PDF appears naturally in turbulence, even without any 
anisotropy. Denote u - longitudinal and v - transverse (vertical) components of the velocity, 
and Ur = u{x + r) — u{x), Vr = v{x + r) — v{x), the velocity increments. Then, we have, 

{ur) = 0, and {vr) = 0, (1) 

and (in inertial range) 

(ul) = -lev, (2) 
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the so-called 4/5-Kolmogorov law [18]. Besides 



(3) 



6 dr 

[lol ]. {20 1. The fact that the first moments vanish, see ([I]), whereas the two third moments 
and (j2]) do not clearly indicate that both the PDF for Ur and the joint PDF for Ur and 
Vr are asymmetric. 

As the Kolmogorov law is derived assuming isotropic turbulence, this asymmetry in 
principle exists without anisotropy. Besides, the scaling defined by the Kolmogorov law 
([2]) does not have any intermittency corrections. Nevertheless, as mentioned above, the 
asymmetry (even without aniso trop y) may be related to the intermittency of turbulence, 
as suggested by [171], [2l| and 



22| |. Still, this connection does not seem to present the 



whole picture. Presumably, as in sheared flows, the anisotropy should also be taken into 
considerations. As mentioned above, even very high-Reynolds-number turbulence manifests 
anisotropy in small scales, which is higher than predicted from dimensional arguments by 
231]. This suggests that there is additional to Kolmogorov cascade transfer of energy from 
large scales directly to small, leading to anomalous events like intermittency and anisotropy. 

In this paper we focus on these events. In particular, we construct experimental joint 
PDF for Ur and Vr to see directly if anisotropy and asymmetry is present in rare violent 
events (responsible for the intermittency). Note that constructing the joint 2D PDF's proved 
to be useful in turbulence, see for example 2J], 25|], 26|]. Even more information about 
the connection between anisotropy, asymmetry and intermittency we obtain from several 
conditional and cumulative averages described in the next section. 



PROBLEM DESCRIPTION 



We will work with dimensionless variables, 

(^2^1/2' r ^^2^1/2' 

and construct experimental joint PDF, p{u'^,v'^) to study both asymmetry and anisotropy. 
This 2D distribution is useful to compare with 2D Gaussian anisotropic distribution, pci see 
Appendix, (USD, (HI. 

The joint 2D PDF gives general information about the distribution. More detailed infor- 
mation which is easier to analyze are provided by different ID distributions and cumulative 
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moments. In particular, we are interested in studding the third mixed moment, ([3]), 

{UrVl) = / UrVlp{Ur, Vr)dUrdVr, (4) 



that is, it is important to find out what part of the distribution contributes most into this 
moment. It was shown before, see j27], that the tail parts of ID distributions satisfactory 
recover the moment. 

Additional information is given by conditional average. 



p{Ur) p{Ur) ' 



or, in dimensionless variables. 



{v'lWr) = [piv'rWrWldv'r = . (6) 

In these expressions we introduced 

^{Ur) = j vlp{Ur,Vr)dVr, ^{u r) = J v'lp{Ur,v'r)dv'r (7) 

Therefore, 

{UrVl) = j Ur{vl\Ur)p{Ur)dUr = j Ur^{Ur)dUr, (8) 

(u'rV'l) = ka = j U r{v'l\u' r)p{u r)du' r = J U r^{u' r)du' r (9) 

The function ^{ur) deserves special attention. It follows from ([7j) that 

j ^{Ur)dUr = {vD, (10) 

that is, in a way, the function is a f ^^-distribution versus Ur- The first moment of this distri- 
bution should not vanish, as it follows from ([8]). This means that the f ^-distribution should 
be be asymmetric. For Gaussian distribution, the function $g can be easily calculated, see 
Appendix, (IT8D . 

We used x-wire data acquired at Brookhaven National Lab. Distance of probe above the 
ground: 35m; Number of samples: 40960000 per component, that is, for longitudinal {u) and 
transverse {v) components. Sampling frequency: 5 kHz. Mean velocity: 5.15076224 m/s; 
rms u- velocity: 1.81617371 m/s; rms f- velocity: 1.3646025 m/s; Taylor Reynolds number: 
10680 (courtesy of Sreenivasan) . As usual, the data are interpreted using the Taylor's 
hypothesis. 

All throughout the paper we process data for closest two samples, that is, r corresponds 
to the smallest distance between two samples. 
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FIG. 1: Experimental joint PDF, core part. The indicated levels correspond to e~^, e~^, e~^, 

of the maximum of the PDF, inside out correspondingly. The contours a) depict experimental PDF, 

and b) - Gaussian anisotropic (fTUj) . 



JOINT PDF 

Figures [1] and [2] present this PDF in different ranges. We (loosely) define these regions as 
core part, Fig. [T], and tail part. Fig. [21 Clearly, Fig. [T] corresponds to the main events, while 
Fig. [2]- to the rare and violent events, - as indicated by the levels given in the captures. 

It is clear that the main events. Fig. [H are not much different from Gaussian. Both PDF's 
are noticeably anisotropic, that is the levels are roughly ellipses with big axises inclined at 
some angles to the x-axis. The positions of the levels are roughly the same. The only 
difference is asymmetry, the latter being absent in the Gaussian PDF by construction. 

The asymmetry of the experimental PDF is evident from the following observations. The 
left edge of the outer level is at u'^. = —4.1, while the right edge of it is at u'^ = 3.8; the left 
edge of the next level is at u'^ = —2.5, while the right edge is at u'^. = 2.4. 

Note that both anisotropy and asymmetry are expected to be manifested by the main 
events. 

Consider now Fig. O We note first that now there is dramatic difference between the 
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FIG. 2: Experimental joint PDF, tail part. The indicated levels correspond to 
e~^, e~^^, e^^^, e^^^ of the maximum of the PDF, inside out correspondingly. The contours 
a) depict experimental PDF, and b) - Gaussian anisotropic (fTUD- 



experimental contours and Gaussian. Namely, the experimental contours are much further 
away from the core values than the corresponding Gaussian. This feature is however an- 
ticipated, simply corresponding to the presence of tails, that is, to intermittency - in both 
longitudinal and transverse velocity component increments. 

On the other hand, other features of the rare events PDF are more surprising. The 
experimental contours in the Fig. [2](a) look roughly similar to those in Fig. [T]^a), only 
rescaled to different values, and - naturally - more ragged. Indeed, we see anisotropy 
- ellipses with axes inclined roughly the same way in both figures. And, what is more 
important, we notice asymmetry in the rare events PDF. The left edge of the outer level 
corresponds to u'^. = —15.4, while the right edge of it is at u'^. = 13.8; the left edge of the 
next level corresponds to u'^. = —10.7, while the right edge is at = 10.2. 

If we characterize the asymmetry by the ratio of the distance from the left edge of the 
level to zero to the distance from the right edge of the level to zero, we will get for both 
figures [1] and [2] - for the external levels the number 1.1, while for the next levels (again for 
the both figures) the value of 1.05. Thus, approximately, the asymmetry is the same in both 



6 



1 0° — 



~ 1 o- 



1 0-4 — 



1 o- 




FIG. 3: Experimental f'^ distribution, compared with its Gaussian analog, and with p{u\ 



typical and rare events. 



THE v'l VERSUS u'r DISTRIBUTION: FUNCTION ^ 



Global features 



Figure [3] presents experimental v' ^ versus u'r distribution. It is compared with its Gaus- 
sian analog on one hand, and with regular experimental distribution p{u'r) on the other. All 
three distributions in the figure are normalized on unity, i.e., 

^{u' r)du'r = j ^ciu r)du' r = j p{u' r)du' r = 1, 

and therefore their direct comparison makes sense. 

It is clear that f '^-distribution has quite extensive tails. Naturally, the Gaussian distri- 
bution is much lower outside its core part. Moreover, the tail parts of $(u'r) are much above 
corresponding parts of piu'r)- This is because the f '^-distribution is a second moment, see 
definition ([7j), while p{u'r) is a zeroth moment of the same distribution, 

p{u'r) = p{u'r,v'r)dv'r- 
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FIG. 4: Experimental f '^-distribution for \u'r\ > t events versus the probability of these events 
P{t) (solid line), compared with its Gaussian counterpart (dashed line). 



Note that the f '^-distribution is much more dispersed at the tail parts as compared with 
p{u'r) in the same areas. Analogous (huge) dispersion of data is observed in conditional 
average, see below Fig. [71 and is discussed in Subsection . We note here only that this 
dispersion is attributed to the presence of intermittency. 

Contribution of the rare violent events 

Additional information about the tails of the v'l versus u'r distribution can be obtained 
from cumulative moments, 

Ir—t roo 
= <l>{Ur)dUr + / <l>{Ur)dUr, (11) 

that is, the contribution of events with \ur\ >t into the moment (f^), see definition of this 
moment in (fTOl) . These moments will be plotted against 

— t roc 

p{Ur)dUr + / p{Ur)dUr, (12) 
-oo Jt 

the probabilities of these events. 
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Figure m presents such a plot. The experimental distribution is compared with its Gaus- 
sian counterpart. For t = 0, all events are presented, and therefore P{t = 0) = 1 and 
(v'ln = 1. If, on the other hand, t ^ oo, then both distributions go to zero. 

The difference between experimental and Gaussian cases is quite substantial. As an 
example (depicted by dashed-dotted straight lines), we see that (v'l) reaches 10~^ fraction 
of its final value with only 2.8 x 10"^ part of all events. According to our estimate, these 
events correspond to quite violent outbursts with t > 27 A As seen from Fig. [3|, these events 
are almost at the very end of the measured tails. The Gaussian counterpart reaches the 
same value of 10"'^ with 4.7 x 10^^ part of events; both these numbers are comparable for 
the Gaussian distribution. 

On the other hand, the (v'l) of the Gaussian counterpart with the same probability 
2.8 X 10~^ as the experimental reaches only 5.8 x 10~® fraction of its final value, these 
numbers being again more or less comparable. 

Asymmetry 

Note that all three distributions depicted in Fig. [3] are constructed with the same reso- 
lution, that is, the bin-size for all three distribution was the same and equal to 0.1. This 
makes it possible to compare the dispersion of the data for different distributions. On the 
other hand, huge dispersion of the f '^-distribution does not make it possible to compare its 
tails. In order to do this, we have to have a more smooth PDF, and that can be done by 
constructing the distribution with a larger bin-size. Figure [5] presents the f '^-distribution 
constructed with a bin-size = 1. 

Figure [5]^a) compares the positive and negative parts of the distribution directly. It is 
obvious that the negative tail prevails over the positive for rather big u'r- That is, the rare 
stormy events are definitely asymmetric. 

The f '^-distribution can be considered as a PDF centered at u'r = ka, see ([9]). That is, 
the distribution is definitely asymmetric, and theoretically, this asymmetry should be present 
at the tails as well. This actually can be seen from the shifted Gaussian distribution $, see 
definition (fT8|) and (122|) . and from the Fig. [5](a), that the negative part of the distribution 
is above the positive for any u'r- 

In order to see if the f '^-distribution is just a shifted distribution, or not, we make direct 
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FIG. 5: Experimental ^'^-distribution. Direct comparison of the left and right wings. 



comparison of its part where u'r — ka is positive with the part where u'r — ka is negative, 
seeO^b). The two parts of the corresponding Gaussian distribution $, of course, coincide 
(and cannot be distinguished in the plot), while the experimental distribution still shows 
asymmetry: The centered negative part still prevails. This confirms the above conclusion 
that the rare violent events are asymmetric. 

Another test for asymmetry is measuring the contribution of rare events, as in Subsection 
. Namely, we will consider cumulative moments, 

(u^f^) I = / Ur^{Ur)dUr + / U.r^{Ur)dUr, (13) 

cf. f llip . Or, in dimensionless form, 

ka{t) = {u'rV'li , ,^ = / * u'Mu'r)du\. + H u' Mu' r) du' r , (14) 

The case t = corresponds to all events, so that ka{t = 0) assumes its final value, and 
at t > it approaches it. The case t oo corresponds to ka = 0. Figure [6] illustrates 
ka{t)/ka{0), or actually it shows how the the third moment is formed. 

Obviously, there is a big difference between the experimental ka{t) / ka{0) and its Gaussian 
counterpart, as seen from Fig. [61 As an example (see the dashed-dotted straight lines in the 
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figure), we took the same probability 2.8 x 10^^ as in Fig. HI corresponding to extremely 
violent events. This time, it reaches 0.03th fraction of its final value. In contrast, its Gaussian 
counterpart reaches the same fraction with 0.01 part of events, - these two numbers (0.03 
and 0.01) being comparable (cf. Subsection ). On the other hand, the Gaussian counterpart 
with probability 2.8 x 10^^ reaches the value of 1.8 x 10^^. Again, these two numbers are 
comparable. 

Thus, the contribution of the rare violent events into the experimental moment ([9]) is 
substantial, as opposed to the "regular" situation presumably illustrated by Gaussian dis- 
tribution. Recall that this (odd) moment does not vanish due to asymmetry of the distribu- 
tions, and therefore this substantial contribution of the rare violent events into /c^-moment 
means that these events also possess asymmetry. 
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FIG. 7: Experimental conditional average (u'^|nV), compared with Gaussian defined in (|2ip . and 
with "ideal" correlation, depicted by dashed-dotted line. 



CONDITIONAL AVERAGES 
Global properties 

It is interesting to note that the the experimental conditional average (f'^|u'r) is much 
above the Gaussian defined in (12T!) . see Fig. [71 If these two variables, v! r and v' r would 
be statistically independent, then the conditional average is unity, well below experimental 
average. The anisotropic Gaussian distribution (for which and v' ^ variables are related, 
and therefore so are v! r and f'^), is still much lower than experimental, as seen from the 
figure. 

In another extreme case (as opposed to statistically independent variables) we have an 
"ideal" correlation, v' r = au'r, where a is a constant, which is, due to our normalization, 
equal to unity. Then, simply, 

{V j.\U r) = U J. 

(cf. with Gaussian conditional average with maximal correlation coefficient C = ±1, from 
f l2ip ). As seen from the Fig. [71 this conditional average is even higher than the experimental 
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value: Of course, any connection between two variables is less than " ideal" . 

Another remarkable feature of this average is its gigantic dispersion, cf. Subsection . 
Perhaps, the simplest way to explain it is to consider the two variables u'^ and f V statisti- 
cally independent. As mentioned above, in this case, theoretically {v'1\u'r) = 1- However, 
experimental measurements would give dispersed values, the statistics being defined by the 
number of events. To be more specific, for small and moderate values of \u'r\ with huge 
number of events, the data of (f^) would be quite close to unity. For big values of \u'r\, with 
only few events the data would be strongly dispersed around unity. In particular, if there is 
only one event in some bins, then the values of (f^) would coincide with v"^. themselves, and 
they would be quite dispersed if the process is intermittent. Qualitatively, the experimental 
conditional average in the Fig. [7] looks as described above, that is, it is smooth for small 
and moderate values of Im'^I, and strongly dispersed for big values. Returning to v'1 vs 
u'r distribution, i.e., function $('uV), "we recall that ^{u'r) = {v'l\u'r)p{u'r), see ([6]), where, 
p{u'r) is not very dispersed function, as seen from Fig. [3l Therefore, the dispersion of the v'l 
vs u'r distribution is explained analogously to the dispersion of conditional average {v'l\u'r)- 

Asymmetry 

The data in Fig. [7] are presented with relatively high resolution, the bin-size being equal 
to 0.1. This makes it possible to see the dispersion and to interpret it as in Subsection . On 
the other hand, it is difficult to study the asymmetry with these dispersed data. Another 
processing of data with bin-size equal to 1 are presented in Fig. [U where we compare the 
negative and positive parts of this conditional average directly. 

We note here that Gaussian expression for conditional average fl2T|) calculated from the 
joint 2D distribution fll6p and used in Fig. [S] is symmetric (recall that it is not known how 
to include asymmetry into Gaussian or some other simple test 2D distribution). For this 
reason, the difference between the right and left Gaussian wings in both figures (a) and (b) 
is spurious. Still, putting the Gaussian conditional average into these two plots seems to be 
useful in order to compare the experimental data with something "regular". 

Both quantities, {v'l) vs u'r distribution. Fig. [5l and {v'l\u'r), Fig. [HI are related. For 
that reason in the latter figure we also check the symmetry in respect to the u'r = ka point 
(in[H](b), analogously to[S](b)), in addition to the symmetry in respect to u'r = point, given 
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in [n](a) andlH^a). In other words, we check if the conditional average {v'1\u'r) is simply a 
shifted (but still symmetric) distribution, or not. 

Close examination of Fig. [8] suggests that there is no substantial asymmetry as a sys- 
tematic trend in this conditional average. Although this is an unexpected conclusion, there 
might be an explanation as follows. As can be seen from Fig. [5l the {v'l) distribution, 
i.e., the function ^{u'r), is asymmetric, and so is the u'r distribution, p{u'r), see |22|. This 



means that, in particular, the negative wings of both quantities are elevated above the right- 
hand wings. According to ([HD, the conditional average {v'l\u'r) is defined as a ratio of these 
two quantities. This suggests that increased value of $('u'r) in numerator is balanced by 
increased value of p{u'r) in denominator, decreasing the asymmetry of {v'1.\u'r)- In other 
words, the latter quantity appears to be less sensitive to the asymmetry than the former 
two. 
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DISCUSSION 



The asymmetry of turbulence has been studied for a long time. Its relation to the 



vorticity production was pointed out by 



281]. The asymmetry of M,,-distribution resulting in 



the 4/5-Kolmogorov law was interpreted with the help of the ramp-model, see [17|], [2l|: As 
(wr) = 0, any compression, with Ur < 0, is as efficient as expansion, with Ur > 0. However, 
the compression appears stronger but rarer than decompression, the latter being weaker 
and longer. Therefore, (u^) ^ 0. This model immediately suggests that this asymmetry 
is related to the intermittency: The compressed rare but strong events are supposed to be 
intermittent. 

The ramp-model is only heuristic, however. It was shown by {29!, [s^, jsil that Burgers 
vortex, embedded into a converging motion, acquires negative skewness, this picture con- 
taining both asymmetry and intermittency. As the ramp-model does not exactly correspond 
;o this picture, it had to be modified. This modification was called the bump-model, see 
271]. Actually, "the bump" corresponds to the Burgers vortex. This model is now two- 
dimensional, - to include transverse velocity increments Vj. into the picture. In particular, 
it automatically explains why the mixed third moment {urv'^), ([3]), does not vanish. 

Fully two-dimensional PDF describing all three items, namely asymmetry, anisotropy and 
intermittency, was discussed in this paper. The PDF, and also some additional conditional 
and cumulative averages seem to support the idea about the connection between these three 
features. In the whole picture, there is additional (to the Kolmogorov cascade) energy 
transport from the large non-universal scales directly to the small scales. We note that 
the Burgers vortex, with a small radius, is generated by a relatively large-scale motion, thus 
making it possible to transfer the energy directly from large scales to small. In the framework 
of the bump-model, it is "the bump" which is generated analogously to the Burgers vortex. 
The scale of the bump is definitely much smaller than the scale of the generating it motion 
(which can be seen from Fig. 5(c) by 23]), and that could explain the direct interaction 
between the non-universal large scales and small scales. 

I thank K. R. Sreenivasan and S. Kurien for the data and for useful discussions. 



15 



Appendix 

The Gaussian joint PDF lias tlie form, 

1 



PG[Ur,Vr 



: exp 



(ur) = (vj.) = 0. Or, for dimensionless variables. 



U;{v;) - 2UrVr{UrVr) + V;{ul) 

2{{ul){v^)-{urV,Y) 



(15) 



1 



u . 



where 



27rVl - C2 

C = {u'rV'r) 



exp 



2u'rV'rC + V'[ 



2(1 -C^ 



(UrVr) 



(16) 



(17) 



1/2 (^2) 1/2' 

the correlation coefficient. Note that in isotropic turbulence C vanishes, and therefore 
the distribution (fTSl) . or (fT6l) . corresponds to anisotropic Gaussian process. Note, how- 
ever, that this distribution is not asymmetric. Indeed, it is easy to show that {u^)g = 

I ulpG{Ur,Vr)dUrdVr = and {UrVDc = J UrV'^PciUr, Vr)dUrdVr = 0. 

For Gaussian distribution (fT5|) . 



$G(«r) 



/2tt 



(^/2)3/2 



(^2)5/2 



Or, for dimensionless variables, we have. 



1 



exp 



/2 ■ 

u' 



2{u^) 



(18) 



271 ^ 2 



(19) 



For Gaussian distribution, according to ([5]) and f|T8l) . we have 



or 



{vl\Ur)G = (Vr) 



{v'1\u'^)g 



1-C' + 



UZ 



(20) 



(21) 



These Gaussian expressions are useful to check in two limiting cases. If C = 0, the two 
variables Ur and Vr become statistically independent, and, naturally, {vl\ur) = {v^). In 
another limiting case of perfect correlation, Vr = aur, and therefore C = ±1, then (f^|ur) = 



2,, 2 



a u 



{vl)ul/{u1). These two extreme cases are indeed con&med by (120!) . 



As mentioned above, the Gaussian distribution (fTSll is not asymmetric. It is not known 
how to construct a simple test 2D PDF containing both anisotropy and asymmetry. We 
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therefore are going to use (fT5l) directly comparing it with the experimental joint 2D PDF. 
However, in all expressions containing ID PDF, instead of symmetric ID PDF resulting from 
integration of the PDF (1161) over dv'j., we will use a simple asymmetric PDF, 7^, constructed 



as a sum of two Gaussian distributions, see 



22| and |27j, so that 



J Irdu'^ = 1, J u'^Irdu'^ = {u'^) = 0, j u'^ Irdu'^ = {u'^) = I, and j u'^^Irdu'^= {u'^). 

Being asymmetric it is otherwise almost indistinguishable from Gaussian ID PDF following 
from (fT6l) . at least visually. 

The Gaussian ^g{x) should be "corrected" as well to incorporate asymmetry. We will 
therefore use 

^g{x')=<^g{x (22) 
instead of $g(^) defined in (fTHl) . Then, of course, 

{u'rV'l)G = J u'r^Giu'r)dUr = ka- 

Note that the difference between the plots of ^g and ^g is only cosmetic. 
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